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ABSTRACT 


A computer program was developed to present shock front properties in the form of isovalues a 
as functions of altitude of hurst and range, This program consisted of a Sachs’ scaling routine and 3 
a log-log interpolation routine, It was found that overpressure, density, and Rho-U decrease with # 
increasing altitude; particle velocity, density ratio, pressure ratio, and shack strength increase 
with increasing altitude; thermal f-ux density and dynamic pressure remain relatively unchanged ; 


with increasing altitude. 
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INTRODUCTION s 
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The standard scientific presentation of weapons effects information is in the form of universal 


; curves a8 in NOLTR-69-88! and EM-%. 2 For engincering purposcs a presentation in the form of 
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; nificant insight into the altitude dependence of isovalue ranges, The obtaining o 
isovalues is a standard and somewhat tedious problem in cross plotting and interpolation. 


Es sealing and interpolation compute= routines were obtained from NOL’ and modified to inc+rporate 
KE logarithmic interpolation. 
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SECTION 2 
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3 DISCUSSION 

2-1, GENERAL, Damage to aeronautical and missile systems within the sensible atmosphere - ; 4 
=< (0 to 100,000 fect altitude) from nuclear weapons effects is expressible in terms of damage distunces & 
ee ¥. 
i for damage environments, such as overpressure, dynamic pressure, thermal radiation, nevtron and & 
Re ee) 
RS X-ray flux, ete. For some systems the dominant or largest damage distances are associsied wi ; 2 
E hydr dynamic shock front properties or thermal! radiation, ‘hese sysiems are generally acre. es 
a aynamic vehicles, such as aircraft, cruise inissiles, and iow-alisiude-point defense missiles, such S 
q as the 3T' s, Nike, or Sprint. = 
a 2-2, SHOCK FRONT PROPERTIES, The hydrodynamic shock front pruperties of interest to aero~ ee 
dynamic vehicles can be defined by a system consisting of a sea-level, 1-kt. pressure-diatance curve 3 
 - and a means of obtaining altitude and yield corrections to this basic curve, s 2 
a 2-3, THERMAL FLUX DENSITY. ‘The thermal flux density (Q) is dofined by the foxzmuia s 3 
9 Q= fWe/R*, where R is the range, f is the thermal partition function, W is the yield, and c isa j x 
Ee constant containing-dimension factors and the transparency factor (T) of the atmogphere (in our . g 
4 case T= 1, perfect transparency). SS 
a . 244. STANDARD 1-KT, NUCLEAR FREE AIR PRESSURE-DISTANCE CURVE FOR SEA-LEVEL t re 
E CONDITIONS, The data for the 1 kt., sea-level cverpressure curve-was provided by two sources: & 
eB above 60 p.s.i., USS9 (DASA 1206) a: 3 

a below 60 p.s.i., NOLTR 69-88. ae 
Be This data is listed as Subroutine RP1271 and plotted in [igure 2-1. 3 
i 2-5, STANDARD 1-KT. NUCLEAR THERMAL FLUX-~DISTANCE CURVE FOR SEA-LEVEL 3 
= CONDITIONS, The thermal flux density formula used was Q = 7.96 fW/R", where 7. 96 includes @ 
Ee =} 
Ss the transparency factor (1), the dimensionless factor a and the unit factor to give Q in terms of z 
z = 
‘a cal, /em.*, The sea-level, 1-kt. thermal flux density is found plotted in figure 2-2. R 
a 2-6, ALTITUDE CORRECTIONS—SACHS’ SCALING, The altitude corrections applicd are the usual : i 
ore ye & 
3g Sachs’ scaling relationships: 3 
Bs AP, _ AP, E 
4 P,P Ei 
3 - When: ze 
a oa 
EE RTP, ] 1/3 _ PP] 12 =: 
‘B W, a ct W, 4 
ag Where: 2 
Fs J = a 
Se AP, and AP, arc ovespressures for point 1 and point 2, respectively, P,; and P. are the respective FE: 
atmospheric pressures, R,; and R, are the respective ranges, and W, and W, are the respective yields. 2 
E 2 F 
E 3 
x 2 
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& s 
. 2-6,1 Subroutine ARDC proviries pressure, temperature, ard density of the standard atmosphere. ; 
a The data was provided by R, A. Minzer, K. S. W. Champion, and H. L. Pond in The ARDC Model 3 
‘ Atmosphere, * 1959, Air Forces Surveys in Geophysics No, 115 (AFCRC-TR-59-267), Air Force 3 
E Cambridge Research Center, August 1959, 4 
4 2-6.2 Using the cverpressure-range curye of Subroutine RP1271 as the kase curve, and the atmos- 3 
ke pheric pressures generated by subroutine AKNC, the ovcrpressures for any altitude and any cor- 3 
2 rected yield (in this exercise W, = W,) can be calculated. 


7 2-7. ENERGY PARTITION FOR BLAST AND THERMAL AS A FUNCTION CF ALTITUDE. The 


yield corrections which are applied are obtained from the more formally corvect radiation- 


x, 


sased inane aldachetiesanncoum dA nae Sentdateernbee 


ito Aa Sateen anSialdateX Katte 


hydrodynamic computer code calculations. These corrections are, unfortunately, functions of yield, 
istance, and altitude; but, for the range of hydrodynamic front and thermal radiation variabies of 
interest in this memorandum, the yicid correction is assumed te be 2 simple function of altitude 


x ealy, The simplifications used in this memorandum are straight-line approximations to the curves 

Bs presented in EM-1, and are applied firs? to the overpressure dependent qualities, and then indeptnd- 

a ently to the thermal radiation. 2 i 
Z 2-7.4 The blast yield reduction with altitude function was taken from the effects manual, *tCapa- 


bilities of Nucleax Weapons, * J January 1968, and the thermal partition function was taken from 
chapter 3, Thermal Radiaticn Phenomena, "* KN-G8-504{R), 26 May 1969. Tuese two functions are 
found pictted as figures 2-3 and 2-4, Also plotted are the straight-line approximations which were 


3 used in the computer program, : 
a 2-7.2 These straight-line approximaticns were used mainly for the purpose of facilitating computer : 
a time, There is no simple functional dependence between altitude and the quartities, and so tu avcid = 
; - feeding in al) the dzia points for every possthie altitude, the straight-line approximation was adopted, Z 
a 2-7.3 For the blast yield reduction curve, the greatest difference occurs around 80 tat,, where 4 
£ there ic s difference of 0.02 (approximately 2-percent error) between the curve and the approxima~ 3 
3 tion. The curve and the approximation agree quite well over the whole altitude range, the erro? 4 
E being less‘than 1 percent over most of the ranze. 4 
s : 2-7.4 The blast yield reduction curve can be compared with a plot of the maximum computed F 
a effective blast yield in DASA 1290, This rict is show in figyre 2-5, The straight-line approxima- 4 
a tion used in the program was fitted ta the curve, not to the envelope of figure 2-5, because in the 3 
E . altitude range of interest there are anly two computations &t 165 Ht,}, and the curve intersects i 
pe these points while the envelope docs not. The function used is cowsidered conservative fer defense Fs 
purposes, as in sure-survival studies cn aircraft, a 
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Figure 2-4, -- The1mal partition function for 1 kt 
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2 ' 2-7.5 For the thermal partition function, the greatest difierence occurs around 20 kft., where 
sy i - there is a difference of 0,01 (approximately 2-percent error) between the curve and the approxima- 
Ble ; 
BS tion, For the majority of altitudes, the error ‘s less than 1 percent. 
g j . 2-8, PRESENTATION OF SHOCK FRONT PROPERTIES (RANKINE-HUGONIOT RELATIONS) 
aaa VERSUS ALTITUL”' CURVES AND THERMAL FLUX DENSITY VERSUS ALTITUDE CURVES. 
iE - i The Rankine-Hugoniot relations‘ enable the calculation of density, density ratio. particle velocity, 
- t pressure ratio, shock strength, and dynamic pressure behind the shock front if the peak over- ° 


ae pressure is given. 
: i 2-8.1 The program presented is not the most efficient program because a slight rearrangement of 
: cards ‘s necessary when changing from interpolation of one quantity to another (i.e. , overpressure 
to dynarcic pressure, deasity to thermal flux density, etc.). 
and straighiforward and preferable to one gigantic, complex program that would just be nine pro- 
: grams added together. This simpie card rearrangement can be seen in the sample programs pre- 
: : sented in appendix II, 
* 2-8.2 The peak overpressure isovalues are presected in figure 2-6, the dynamic pressure iso- 
i . values in figure 2-7, the particle velocity in figure 2-8, the density ratio in figure 2-9, the parti- 
r cle velocity x density (Rho-U) in figure 2-10, the pressure ratio in figure 2-11, the shock strength 
in figure 2~12, and the density in figure 2-13, 
2-8.3 Ali quantities were plotted from sea level to 109-kft, altitude, except for density, which is 
> plotted only from sea level to 40 kit. in the most advantageous cases. The reason for this can be 

: discerned by doing a few simple calculations. 
: 2-84 The formuls for density is DENSITY = 0.0764575 x DRATIO x ean where POVPA= 
. : SSuae* Since DRATIO = 1 at sea level, the minimum end maximum values of DENSITY are 
ed entirely dependent on the theoretical limits of OVP. 
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If OVP « 0, DENSITY 
If OVP + ©, DENSITY 
At 50 kft;, DRATIO 


Therefore, if OVP - ©, DENSITY = 0.076475 x 


= 0, 076475 


= 0.076475 x6 = 0.45885 
0. 012709 
0, 075675 * 


0, 014799 
9, OTEATS 


x 6x 0, 670254. 


2-8,.5 Since the highest possible value of DENSITY at 50 Kft, is less than the lowest possible value 
of DENSITY at sea level, it is tmpossible ta produce a plot that will have an altitude dependence 


from sea leve? to greater than 40 kit, 


2-8,.6 The thermal flux density isovalues are presented in figure 2-14. 


‘ Glasstone, S. (ed.). "The Effects of Nuclear Weapons." 1964, 
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Figure 2-6, ... Free-field overpressure (p, s, i.) a8 a function 
of altitude and range, 
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Figure 2-9, -- Density ratio as a function of altitude and range, 
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. Figure 2-12, --Shock strength as a function of altitude and range. 
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2-9, YIELD CORRZCTIONS, Distance upproximations for values of yield Gther than those presented 

in this memorandum may be obtained by muitiplying the range by the cube root of the yield for hydro- ° 
dyramic properties and the square root of the yield for thermal radiation. Discrepancies in the 
yield-correction factors are masked somewhat by the cube root and square root denendence, but 
very lerge yield variations may require more accurate estimates of the yield-corrections {actor as- 


prescoted in EM-~1, 
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Figure 2-14, -- Thermal flux density (cal, /em.?) as a furction of 
altitude and range. 
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SECTION 3 
CONCLUSIONS 


3-1, It is found that overpressure, density, and Rho-U decrease with increasing altitude; particle 
velocity, pressure ratio, shock strength, and density ratio increase with increasing altitude; ther- 
mal flux density and dynamic pressure remain relatively unchanged with increasing altitude. 

3-2, This typ of presentation should be very useful in determining approximate sure-kili and/oz 
sure-survival ranges for the various zerodynamic vehicles, By placing the various isovatues on a 
single chart, a sure-kill envelope could be drawn for the range 0 to 100,000 feet. For exemple, if 
the sure-kill criteria is an overpressure of 2 ¢,3,i,, a particle velocity of 100 ft. /sec., or a ther- 
mal flux density of 20 cal, /cm,?, and a sure-eurvival criteria is an overpressure of 1 .8.3., & 
particle velocity of 50 ft. /sec., or a thermal flux density of 10 cal. /c.n,”, then the determining 
sure-kili/surs-survival range for a 1-kt, explosion is entirely dependent on the particle velezity, as 
shown in figure 3~1, 
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Figure 3-1. --Free-ficld overpressure (0. 8. i.) thermal flux donsity 
(cal. /om,?}, and particle velocity (ft. /sec.) for 1 kt, 
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. . APPENDIX I 
Z INTERPOLATION SCHE ME 
a . 
e i The interpolation scheme used is a straight-line approximation scheme; i.e,, it is assumed that 
L- ‘ in a well-behaved region of a curve, two points can be chosen, and the straight line defined by theze & 
F : points can be used to generate a value of R that will approach the desired value of R, By repested : 
EF use of this scheme, the desired value uf R can be approached to within any desired accuracy. The % 
Fe & 
He program actually used an accuracy of 0,1 percent in R. BS 
fe An initial value of R, is chosen, and then R, is chosen to be equal to 2R;, The scheme is as ‘ 
Fs ‘ 2 
follows: 4 
BE 3 
slope = Ry) -f ¥ 
E. Ry - R 3 
es : s 13 
—— ‘ obviously, slope = fR) {(Raesired) 3 
& ; Ry ~R Z 
1 new & 
P - aa 
therefore, R= Rt A f(Facsired 7 ACh) 3 
j 3lope R 
ioe If | Row Ri| = 0.001R,, then Rew is indeed R destred* “ {R. oS Rif > 0,901R,, then 4 
: R, Rr R 4 
me “te, Raa R 
new 


and another value of R is generated until the condition of [Row - Rif = 9, 0OLR, is met, 


. This method of interpolation is valid as long as the foliowing conditions are met: 
Q) £(R) = f @,) for R, = R, 
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(2). £(R) is continuous over the specified range 
‘ ’ (3) £7 (R) is continuous over the specified range. 
: Obviously, the clorer £ (R) is to a straight line, the faster this acheme will work. For this 
ES : Zeason, since the base curve overpressure plot, which is the main curve for calculating most of cur 
E cae quantities, has a dependence on R similar to Pia acai and Q<R™, & log-log plot was used; that 
EB: is, instead of slope = ae , the formula slope = | 
. i ah as s ee was used, ; 
e The use of this latter formula is more than just a converience, it is a necessity when R % 
Ee ‘ approaches 0, At this pcint (R= 9), the necessary conditions that f(R) and £“(R) be continuous no 4 
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‘onger hold. By "flattening" the curve with the use of log £(R), the discontinuity at R = 0 can be 


much more closely approached. 
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